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Abstract 

In this paper we continue the investigation of partition functions of critical systems 
(~| , on a rectangle initiated in [R. Bondesan et al, Nucl.Phys.B862:553-575,2012]. Here 

we develop a general formalism of rectangle boundary states using conformal field 
theory, adapted to describe geometries supporting different boundary conditions. 
We discuss the computation of rectangular amplitudes and their modular properties, 
presenting explicit results for the case of free theories. In a second part of the paper 
we focus on applications to loop models, discussing in details lattice discretizations 
using both numerical and analytical calculations. These results allow to interpret 
geometrically conformal blocks, and as an application we derive new probability 
formulas for self- avoiding walks. 

1 Introduction 

The study of conformal field theories (CFT) on a rectangle, while relying on a well 
known theoretical framework, offers interesting technical challenges, and had not been 
considered much in the literature until quite recently. In the last few years, the topic has 
however received increased attention, both in the context of open string field theory[22, 
23,37], and in the context of condensed matter physics applications. The latter include 
study of quantum quenches for systems with open boundary conditions [6, 15], or the 
study of transport properties in disordered network models [4,10]. 

For a CFT on a rectangle, the topology is trivial, while there are sharp corners, and 
in general four different boundary conditions on the four edges, see figure 1. 

The most natural way to calculating properties — such as partition functions — in this 
geometry relies on a concept of "fully open" boundary state, which was introduced in 
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Figure 1: The partition function for a CFT on a rectangle. Different boundary conditions 
are imposed on each side, and boundary condition changing operators are inserted at 
the corners Zi. 

[3]. In the special case where all boundary conditions are the same, this state was 
studied in details in [3], and related to remarkable algebraic properties of the Virasoro 
algebra. In the same reference, the correspondence between the CFT and lattice model 
descriptions of the rectangle geometry was also considered in details. The present paper 
is a continuation of this work, where we now consider the case of different boundary 
conditions. 

The first two sections delineate the general formalism. In section 2 we introduce the 
most general fully open boundary states \B^^) and discuss their general expression in 
terms of basis states which are the natural analog of the Ishibashi states [24] in 
this context. In section 3, we study amplitudes (partition functions), and compare the 
approaches using the open boundary states with the approach using conformal mappings 
and conformal blocks. Various checks and related results are considered in section 4 
where we discuss free theories. 

We then turn in section 5 to our main motivation for this work, which is the calcu- 
lation of partition functions for loop models of statistical mechanics with lines inserted 
at the corners, for an example see figure 2. These objects have potential relations with 
transport properties of network models as well as potential probabilistic applications. 
They also provide a natural tool to investigate degeneracy issues when the correspond- 
ing conformal field theory becomes logarithmic [32,33]. Although the conformal field 
theoretic description of loop models is not fully understood, we are able to find out the 
detailed connection between the general considerations of the first part of this paper and 
boundary conditions defined geometrically. This allows us to interpret conformal blocks 
in terms of loops, and calculate the desired partition functions. We find in particular 
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Figure 2: A schematic illustration of the partition function for a loop model where we 
insert i lines in the bottom left corner and j in the bottom right one, and we impose 
that s = i + j lines propagate vertically. 

that coefficients in geometrical amplitudes for percolation, dilute or dense polymers are 
then related with indecomposability parameters. 

A few final remarks are gathered in the conclusion, while some technical points are 
addressed in the appendices. 

2 Solutions of the gluing condition for semi-rectangular 
geometries 

A boundary condition in conformal field theory has to preserve conformal symmetry. 
This is encoded in the so-called gluing conditions for the stress tensor of the CFT, 
gluing left- and right-moving modes. In the familiar case of a theory defined on a disk, 
this implies that 



where \BP) is the boundary state describing the boundary of the disk. \B^) is an element 
of the space of states F of the bulk theory which we assume to decompose onto Virasoro 
irreducibles Vi with multiplicities Nij, F = NijVi iS" Vj. The condition (1) can be 
solved in each direct summand separately and the unique solution in (g) is given by 
so-called Ishibashi states [24]: 



with \i;N) {\i;N)) an orthonormal base of Vi (Vi). 

Using the transformation law of the stress tensor, also more complicated boundaries 
can be described by gluing conditions. In [3] the geometry of the bottom of the rectan- 
gle, see figure 3, was considered. In this case not only left and right movers are glued 
together, but also positive and negative modes within a chiral algebra. Further, confor- 
mal anomalies coming from singularities of the stress tensor at the corners will show up 



(L„ - L_„) \BP) = 



(1) 




(2) 
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Figure 3: The semi-rectangular geometry defining the boundary state. 

in the gluing condition. Consider the possibility of having different boundary conditions 
whose change is mediated by fields inserted in the left and right corners of the bottom 
of the rectangle, of weights respectively hi and hr- Then, denoting the boundary state 
\B°), the gluing condition reads: 

(l„ - L_„ - 2n (hi + i-rfir)) IB") = . (3) 



where 



h = 2h-^, (4) 



c being the central charge. 

The rectangle bottom in the (^-plane, with corners say in = —1, 0, can be mapped 
to the region V = {z G C : Q{z) > 0,\z\ > 1}, left of figure 4, by z = e-'""^ . This 
geometry is the appropriate one to discuss boundary states in radial quantization, and 
clearly the gluing condition (and thus the boundary state) is not changed with respect 
to that for the bottom of the rectangle. In [3] the operator implementing algebraically 
the mapping from the upper half plane to the geometry with corners T>, mapping z = ±1 
to w = ±2, see figure 4, was found. It has the following product formula: 

Gv = lim e"W=T^-2^ . . . e-i^-^e'^-^ . (5) 

N^oo 

It follows that when hi = hr = 0, homogeneous boundary conditions, the boundary 
state is: 

IB") = Gv\0) (6) 
= |0) - L„2|0) - ^L_4|0) + + .... (7) 

We will comment about the normalization of the state \B") when we compute the am- 
plitudes involving boundary states below. In the following we will continue to use the 
normalization fixed in equation (6). 
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Figure 4: Mapping of semicircular region D to the upper half plane H. / is the Joukowsky 
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In this paper we focus on the case of different boundary conditions on the sides of 
the semi-rectangle. To make things concrete, call the condition on the left side a, that 
on the bottom b and that on the right c. Boundary condition changing (BCC) operators 



-which can be identified with the lowest eigenstate of the Hamiltonian on 



a strip with boundary conditions a and b on the sides — will then sit at the left and 
right corners. We denote the boundary state we want to determine \B^^), which is thus 
associated graphically to: 



(8) 



UX = ^zi,zi] 
is defined by 



(j){x) • • • is a chain of arbitrary bulk and boundary operators, |-B^c) 



{X)v = {o\x\b: 



^ ) 

act ' 



(9) 



where (X)x) is the correlator in the geometry T>. The operator Gx> introduced above by 
definition implements at the operator level the conformal map from T) to H, so that we 
have: 



(X)^ = (0|X,^f (-2)0^(2)10), 

where we defined X = G:^^XGv (if Gv = e^n^r.L„^ ^^len G^^ = e"^" 
right term is a correlation function on the upper half plane. Since the out-vacuum (0| is 
invariant under this mapping, one has 



(10) 



and the 



\B' 



G'i,</.f (-2)</.f(2)|0). 



ill) 



We now recall the boundary OPE of two boundary fields (see e.g. [28], x > y): 



ia,b I 



(12) 
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The boundary OPE coefficients Cf^^ are part of the known data of the CFT [17,34]: 



(13) 



where the F-matrices relate conformal blocks in different bases [31]. We will see explicit 
examples below. Note that from taking different OPEs of the three point function 



l^ab^bc^ai^ One has the cyclic relation 



Cabcr-iaca 
ijk^kkO 



C'bca /^aba 



(14) 



Further, the normalization of the identity </>o (0 is the label for the identity) is fixed by 
demanding 4>(f'4'i'' = (pf^^ which implies 



f^aab 



for every a,i. More properties can be found in [28,34]. 
Eq. (11) can then be written as: 



r6 \ _ \ riabc I r-iaca\ s \ 
^acl — ^ijs V ^ssO \i jl 



sO\i jl ' 



(15) 



(16) 



where we defined the basis states 



by 



k j/ 



n>0 7l-n 



ths—hi—hi 



n>l -yhn 



(17) 
(18) 



where 7 h n is a partition (71,..., 7^) of n, and we used the shorthand L_ 



The coefficients /S^'J appear in the OPE of (bulk) chiral fields 



in the term corresponding to the descendant of s, L-^(f>s (We use the notations of [18], 
section 6.6.3 of this reference, although we simplify the notation by writing 7 instead 

7^gQ in eq. (16) is needed in order to define the basis state as 

without any information about the 



of {7}). The term 

an element of the Verma module of the primary 
boundary conditions. Alternatively one could remove the factor y/C^^Q and define the 
basis state as an element of the Verma module of 4>s'^, where {(ps'^lcps'"") = C^gQ. 

The basis states span the space of boundary states solving eq. (3) and are the 
analogues for the rectangular geometry of Ishibashi states (2). For determining the 
coefficients (5's in eq. (18) one could work out the algebra in (17) or directly solve the 
constraint (3) for different values of n. It is well-known however that the combinatorics 
of the Virasoro algebra is complicated and does not allow to get a closed formula for 
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the coefficients in the above expansions, be it fS^'J or 5f'J. An exception is the case 



of free theories, when the Virasoro generators admit an expression in terms of bosonic 
or fermionic oscillators, see the discussion in section 4. In general, there is therefore 
the following difference between Ishibashi states for the cylinder geometry and the basis 
states for the rectangle: rectangle basis states depend on the boundary operators inserted 
at the corners, so on their fusion, which does not give access to an explicit expression. 
Instead Ishibashi states are known explicitly from eq. (2). 

The state \B^^) is directly given by a single basis state if the fusion (j)^'^ produces 
only one channel. This is of course the case of homogeneous boundary conditions, when 
no BCC operators are present. 



where \(j)Q) = |0), and we used (15). This is also the case when there is only a single 
BCC operator. Indeed, if o = 6, calling the BCC operator (j)^'^, one has 



since -L_i generates translations and \(j)s) inserts the operator at the origin. C^^g "^^^ 
fixed by (15) and could be set to one, see [34]. 

The discussion presented so far solves implicitly the problem of finding the boundary 
states for the rectangular geometry. In the next section we will compute amplitudes 
involving these states and give more explicit results. Before moving on, let us define the 
conjugate states of eq. (16) 




(19) 




(20) 



./C^GvUW) = 7C^G7?exp(2L_i)|<A,), 





(22) 



s 



which are associated graphically to 



h 




(23) 



(*/| is the conjugate of and we introduce also the following graphical notation for 
the basis states: 




(24) 
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3 Amplitudes and conformal blocks 

We consider a rectangle of length L' and width L in the z-plane, with boundary condi- 
tions a, b, c, d, so that BCC operators (f)'^"'^ (zi) , (j)'^'^ {Z2) , (f)'!^'^ (z^) , (/'f'^ {Z4) sit at the corners, 
see figure 1. The CFT partition function is by definition: 

Z [ate) =^lp^=Zo{T){ct>f{z,),(P';'\z2),4'{z3),4''{z,)) , (25) 

where Zq is the partition function for homogeneous boundary conditions and the corre- 
lator is evaluated on the rectangle. The arrow in the picture is the direction of imaginary 
time in a transfer matrix picture, which will be developed below, and r = iV /L. 
We recall first that Zq has been computed in [26]: 

a 

Zq{t)= «lrl« = L^/\"^/2 (26) 

where r] is the Dedekind eta function. Definitions and properties of the r] function 
are collected in appendix A. In [3], we showed that the following relation between the 
boundary state (6) and this partition function exists: 

q being the relevant combination of L' /L appearing in the strip transfer matrix: 

q = = exp(7rzr) . (28) 
From relation (27) we note that the state \B°) is not normalizable: 

{CXD if c > 0, 
if c < 0, (29) 
1 if c = 0. 

A similar situation happens for the Ishibashi states in the cylinder geometry. However in 
the case of cylinder states, one can fix a possible multiplicative constant in the definition, 
which is not fixed by the gluing condition on the disk (1), by requiring the equality of 
the cylinder amplitude to the annulus partition function, see [8]. We do not know how 
to fix a possible multiplicative constant a in the definition \B°) = aGx> |0). This issue 
is related to the very definition of a partition function on a rectangle, containing the 
factor L^/^ (or L'^ in the general case, see below) in (26), which depends on the unit of 
measure of length we adopt. 

We now go back to the computation of (25). We choose the following fusion channel 
of BCC operators (here depicted on a segment): 
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This corresponds to fuse fields along the bottom and top segments of the rectangle. The 
partition function has the following expression in a transfer matrix formulation in terms 
of boundary states: 



Z { a'^c 



h 



^ ) 
act 



T w(i,j,k,l) \ ^ /^adc/^abc/~iaca li ii'-Lo— c/24| s \ 
^ 2_^'~^ils '~'jks'~^ssO\ slQ \j k) ■ 



(30) 
(31) 



The modular weight w{i,j, k, I) appearing in the term L^^^'^'^'^'^ is given by the anomalies 
at the corners [12]: 

w{i,j,k,l) = ^-2{hi + hj + hi + hk). (32) 

The building blocks of the partition function (31) are amplitudes involving basis 
states: 



(33) 
(34) 



n>0 



In practice using the expression for boundary states (17) one can only obtain the first few 
coefficients c„ of the power series above. However closed expressions can be computed 
using CFT correlators of the BCC operators, since by definition one has the following 
relation: 

= Z,{T){Uzi)H^2)<Pk{zM^,))s , (35) 

where only the channel s in the correlator is considered. For computing a correlator of 
the fields on the rectangle, we map it to the upper half plane w 



h{zi) ■ ■ ■ MZn)) = {Z'{wi))-^^ ■ ■ ■ {z'{Wn)r 



\_{wi) ■ ■ ■ (j)n{Wn)) ■ (36) 



This mapping is given by a Schwarz-Cristoffel transformation: 

'^^ 1 



(37) 
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It maps the points on the real hne (— 1/fc, —1, 1, 1/k) to the corners of a rectangle with 
L = 2K and L' = K', {-K + iK',-K,K,K + iK'), where K is the complete elliptic 
integral of the first kind of modulus k and K' the complementary one, see figure 5. 



-K + iK' K + iK' 





-K K 

Figure 5: The Schwarz-Cristoffel transformation of eq. (37) sending the the upper half 
plane to the rectangle. 

The Jacobian J of this transformation is singular at the corners, but when the 
boundary operators sit at the corners too, eq. (36) continues to hold since as the argument 
of the field approaches the corner in Zc, the field itself goes to zero as {z — Zc)^(j){zc), 
where h is the dimension of the field. Then we look at the Jacobian as the coefficient of 
the first term in the expansion around the position of the corners. The result is 

J — Jj-hi-hi^-y _ f^2-jh^+hj+hk+hi ^gg^ 

Thanks to global conformal invariance, a chiral correlator can be written as [18] 

{(t)i{wi)4>2iw2)(t)3{w3)(p4{w4)) = Yl'^ij' '^(0 ' f^ij = ^ i^hkj - hi - hj , (39) 

i<j \k=l I 

with being the anharmonic ratio (wij = Wi — Wj) 

C = ^^^. (40) 
W13W24 



We now let (twi, ^2, 1^3, W4) = {—1/k, —1, 1, 1/A;) and we find 

JJy^Ai^j ^ ^j^4:(hi+h4)+h2+h3 |'2(i _ ^2^-j-('ii+'*2+/i3+/J4)y/3 ^ ^^^^ 

Expressing everything in terms of r, using eq. (199), we have: 

237rV(2i^)-') G(C). (42) 

We recognize in the above formula the term length L = 2K raised to the correct modular 
weight (32) coming from field insertions. G{C) in general is a linear combination of 
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conformal blocks -^/^.^^(C); corresponding to the fusion channel s, and in the end eq. (33) 



becomes: 



4'^-''=''(r) oc r/(T)-2-(^'^>'')j-?^ (1 _ ^(^)) , (43) 



where the anharmonic ratio as a function of r is 

k-lV (Oiir)^^ 



Let us compute the constant of proportionality by demanding the normalization of the 
amplitude to be cq = ^^'•-'y>^^+^^j+^k+hi) (34)^ according to eq. (18). When g — ^ 0, 

1 — C(t) —7- 0, and using that our conventions are such that 

^u-A^) ~ z^'-^i^^^+^^+^>'+^'^ , (45) 

we have the following small-g expansion of the right hand side of (43): 

Then the exact relation is: 

3.1 Dualities and modular covariance 

Duality relations act in the space of conformal blocks [31]. The F-duality states the 
associativity of fusion product, relating conformal blocks singular near = to those 
singular near C = 1- 

-^4h(C) = E^-[ j i]^i-Jki^-0- (48) 

r 

Another transformation is braiding, exchanging the position of and cpk'- 

■^4h(C) = E^-[;- ■ (49) 

F- and i?-matrices satisfy several identities [31]. The relation between F-matrices and 
boundary OPE coefficients Cfj'^ has been pointed out in eq. (13). 

Let us remark now that the space of conformal blocks furnishes a representation of 
the modular group T = PSL(2, Z) with the following natural action: 



M o t{C{t)) = t{c{Mt)) , Mer. 



(50) 
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Indeed one can relate C,{Mt) to C(r) using the explicit expression of the anharmonic 
ratio in terms of the modular parameter, eq. (44). For the generators 5" : r — )• — 1/r and 
T : r — )■ T + 1, one has 

1 



C(5t) = 1 - C(r) ; aTr) 



so that the action of T and S on conformal blocks is 



r 



j k 

i I 
j k 



C(t)' 



J^i-jkiC{r)) 
-^fc;,/(C(r)). 



(51) 

(52) 
(53) 



This action transfers immediately into modular covariance of rectangular amplitudes 
using eq. (47). Modular inversion S acts in an obvious way on the rectangle by inter- 
changing its length and width. Modular translation T instead deforms it to a rhombus, 
and its physical meaning could be understood by using the action given in the equation 
above as braiding of the boundary fields. Let us see then in more details how T acts on 
the rectangle. S'-action allows to relate amplitudes of rectangles of inverse aspect ratio 
as 



5o4'J>''(r) =^i«'^(-l/r) 



-w{i,j,k,l) 



j k 



where we defined 



with 



^i«>(-l/T) = ('/'||^«-^/24|/A 



exp(— vri/r) 



(54) 
(55) 

(56) 
(57) 



Note the appearance of the modular weight w{i,j,k,l) (not present if considering the 
transformation of the full CFT partition function Z, because of the relation (31)) from 
the transformation of r], see Appendix A. To see how T acts, we consider its action on 
the amplitude of a rectangle of inverse aspect ratio 



To 4'J''=''(-l/r) = 4'J''=''(-1/(t + 1)) 

_ ^i-n-w{i,j,k,l)/3^^ _^ i^-w{i,j,k,l)^^^-^-2w(i,j,k,l) -ps ^ ^ 



(58) 
(59) 



^i-!T'w(i,j,k,l) /3 



r+1 



—2w(i,j,k,l) 



j k 



^;'^'^''(-l/r) . (60) 
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Note that the amphtude 1/r) can be thought of as the r-channel amplitude of 



a rectangle whose bottom, where operators and (/)^'^ sit, has been twisted by tt: 




Recall that for the torus a Dehn twist by 2tt at fixed time is exactly T; here we have q 
so the twist is the half. This should be the way to understand modular translation of 
rectangular amplitudes. In particular if the twisted amplitude is equal to the original 
one, as happens when k = j, we expect to have invariance (in the sense of equation (58)) 
under modular translations. 

More generally modular transformations are powerful tools for constraining the ob- 
jects we want to compute, and their usefulness are well appreciated for the case of a CFT 
on cylinders and tori [7,8]. Let us see what symmetry constraint one has on the rectan- 
gle^ and how the basis state we have introduced solves it. The natural constraint comes 
from imposing the consistency of the partition function upon switching the direction of 
imaginary time from vertical to horizontal: 

d . d 



a] 1 [c — a 



From eq. (31), the 1. h. s. is 

a] T Ic riiifj.i.fcJ) \ ^ /~<adc/~<abc/~iaca Ai,j,k,l/^\ 



(61) 



(62) 



while the r. h. s. is 

d 



j^>w{i,j,k,l) \ ^ ^dab ^bcd^dbd ik 1\~Lq- 
/ J I J s kls ssO \ s iy 



■c/24| 



(63) 

j^lw(i,j,k,l) ^-^{hi+hj+hk+hi) ^dab^bcd^dbd^^_-^^^s^-2w(i,j,k,l) ^g^^ 

s 

X J-^.,,(1-C(-1/t)) 



T w(i,j,k,l) \ ^ \ ^ /~idab /~ibcd/~idbd 77 
^ ^ijs ^kls ^ssQ^s'^ 



(65) 



^ The action of the modular group on a rectangle has been discussed in [13, 27] in the case of 
percolation. 
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Equating the coefficients of the amplitude we get the following consistency equation: 



This equation is simply stating the associativity of the fusion product of boundary fields, 
the so-called sewing constraint coming from crossing symmetry of correlators [28]. It 
has been shown in [34] that its solution is given by (13). Then, since we defined directly 
our boundary states in terms of physical boundary OPEs, consistency under modular 
inversion does not give any additional constraint. Recall that for the cylinder geometry 
a similar argument gives the so-called Cardy states [8] . 

3.2 Examples 

We now show how to construct boundary states and associated amplitudes in some 
explicit examples. The case of homogeneous boundary conditions — no BCC operator 
insertions at the corners — has already been discussed in details in [3], and we here 
concentrate on the change of boundary conditions. 

3.2.1 Two operator insertions 

The simplest case to deal with is that of two different boundary conditions around 
the rectangle, when two BCC operator are inserted in two corners. This computation 
amounts to transforming the two-point function of the BCC operators in the upper hal 
plane to the rectangle geometry with operators in the corners and has already appeared 
in [26] . In this case only one channel can propagate and our formulas of section 3 are of 
course simplified considerably. For definiteness let us consider first the case of boundary 
condition a on the bottom and b on the other sides, and call 0"'^ the BCC operator. The 
partition function is: 




(66) 



s 




b 



a 



(67) 



c/4-4:h^aba yiO,i,i,0 



(68) 
(69) 



Using eq. (47), we express the amplitude using the following conformal block: 




h/3 




8h 



(70) 
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to have: 



_ ^abaA—2hTc/A~Ah—c/48 

— *-^mO 4 q 



\-Sh 



1 + 



32h 



+ 



c(6 + c) 



8(2 + c)h + 128/1^ + . . . 



(71) 
(72) 

(73) 



Note that as expected from the symmetry of interchanging the boundary operators, this 
partition function is invariant (up to a phase) with respect to modular translation. Of 
course it is not invariant under modular inversion S, and indeed under S this partition 
function transforms into that for operators inserted in the left or right corners, which is: 

b 



Z { a f^b 



Z { b a 



b] T a 



T c/i-ih I -Tjb |-Lo-c/24|d6 \ 
^ \i^ab\(l l^ab) 

(^a6a^c/4-4h^(^)-c/2+16/.^(^/2)-8^ 
^aba Tc/4-4:h-c/4S+h/2 



c + 8h(8h-l) 
l + 8hq + ^- ^-q + 



(74) 
(75) 
(76) 
(77) 

(78) 



If instead the operators are inserted in non-adjacent corners, say in the top-left and the 
right-bottom one, we then have 



Z [b^a 



blf^a: 

a 



Z{bta) 



a 

bVV\a 



T c/4—Ah / Tjb i-Lo— c/24|Da \ 

= L'=/4-^^ciro'^°'''°'V) • 

The following conformal block 

^„,,o(i-o = (c(i-c))''' = 2«'»(5e^|5jj^ 

contributes to the amplitude, and we have: 



8h 



Z b 



C^^^L^/^~^''r]{Ty^l^-^^r]{T/2f'^7]{2Tf^' 

ribab T c/4-4/i„-c/48+/i/2 



, . C + 8/l(8/l-l) 

1 - 8/ig H i -q + 



(79) 
(80) 
(81) 



(82) 

(83) 
(84) 



The above results can be also obtained by computing the overlap of the boundary 
states. The amplitude *'*'°(t) involved in the partition function of eq. (67), is defined 
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as in eq. (33), so it will be given by the overlap of the boundary state for homogeneous 
boundary conditions, see eq. (19), and whose first levels read: 



We have checked that {B^i^\q^°~^^'^'^\B^f^) computed using the above expression of bound- 
ary states reproduces the small-g expansion of eq. (73) up to order q^^. The boundary 
states building the amplitudes of eq. (74) and (79) are given by a single basis state as 
in (20). Eq. (74) is given by the overlap of: 



with itself. Instead the partition function in eq. (79) is given by (^j*|^^''~'^''^^|j*o)- (''/I 
is computed from (^^''l by sending 2 — )■ —2 in the explicit formula (20) and this gives 
the minus signs present in equation (84) with respect to (78). Again we have verified 
the agreement of the series expansions using the explicit form of the boundary states 
up to level q^^. Note that the small-(; expansion (73) contains only even powers of q, 
as expected from left-right symmetry of the boundary conditions. Indeed odd powers 
are related to descendants at odd levels changing sign under inversion of the coordinates 
z — )• —z, and thus are not allowed to propagate. 

Furthermore one can check the agreement of these results with computations for free 
theories presented in section 4, eqs. (105), (106) and (107). 

3.2.2 Four operator insertions 

We now turn to the case of four insertions of BCC operators, and we consider the simple 
situation when the operators inserted are all degenerate at level 2, so that the dimension 
of each operator equals hi^2 or /i2,i. We make this choice since an explicit and simple 
expression for the conformal blocks of these fields is available [18] (actually it suffices that 
just one of the four fields is degenerate at level 2 to find a closed form of the conformal 
blocks, but we treat the case of all fields degenerate at level 2 for easiness of notation). 
For fixing ideas we can think of the Q-states Potts model with fixed boundary conditions 
on left and right boundaries and free on the top and bottom sides. Clearly we have to 
distinguish two situations: either the fixed spins are in an equal state a or they are 




(85) 




(86) 
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different, say in states a and a'. Tliese partition functions liave been already discussed 
in details in [11] and here we reproduce the results and comment about boundary states 
and modular properties. Introducing the label i of the Virasoro modules V(i (or 
V(i+i^i)) in Kac notation, one associates to free boundary conditions on the strip the 
label 1, to fixed equal 0, and to fixed different and 2. This is consistent with the 
known fusion rules [11] since 1 (gi 1 = © 2. Then concretely we would like to compute 
the partition functions 

Z{o\a) =0]Tla^ a = 0,2. (87) 

If a = 0, at the corners on the bottom there are BCC operators (j)^'^ and whose 
OPE coefficients are clearly C^'l'^ cx 6sfl so that only the identity channel propagates. 
For a = 2 we insert instead (^^'^ and </)}'^ whose OPE coefficients are C°'J'^ oc 5s,2 so that 
only the tp^-daannel is left. Calling h = hi and Aa = Aa'^'^'^ , we then have 

Z {o\a) =Af^L''/^-^''Aa. (88) 

The normalization A/"" can be written in terms of OPE constants, and its value will 
be fixed later. We turn then our attention to the computation of the conformal blocks 
-^111 1 present in the expression of the amplitudes of eq. (47). This is a standard exercise 
[18]. From the null state condition at level 2, the function Q"', 

g"(C):=(C(l-C))''/'-^f,i;i,i(C), (89) 
satisfies the hyper geometric equation 

C(l - C)Q" + {7 - (1 + a + /3)C} g'-apg = 0. (90) 

with parameters 

a = -4h,f3 = ^-^h,^ = ^-^h. (91) 
The roots of the indicial equation 

ri = , r2 = 1 - 7 , (92) 

label the two solutions of the differential equation (we take the ~ solutions), defining 
the conformal blocks of the identity 4>o, ^'^(C)) and that of the operator (p2 of weight 
/i2 = 8/1/3 + 1/3, ^^(C): 

.»(C,^...(l-|.-..M;c) 



^9,,^ .8^4-1 „ /'I 4/l 4/l 2 8/l 4 ^, , 

g2(^)=^-+3 2Fi (--—,— + -;— + -;C) . (94) 
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The validity of these solutions is for 7 not an integer. 
After a little algebra the partition function reads 

z{o\a)= M''L^i^-^^r]-^i'^ef'^g''{i - C) , 



(95) 



where we reabsorbed all constants in the definition of A/"". The relevant F-matrices for 
the conformal blocks involved can be computed using standard hypergeometric identities 
relating the conformal blocks of the correlator ((/)i(/>j(/)j(/)fc), allowing explicit checks that 
they behave as expected under modular inversion. We report here for future convenience 



the F-matrices Fi 



■ u 



1/2 1/2 
1/2 1/2 



-^00 -^02 
-^20 -^22 



1 



r(4/i+i) 



V 3 +3, 



3 37^3 3 j 

r(|-f)r{-4.) 
_i 

/3 



(96) 



where we defined /3 = 2sin (7r(8/i + l)/6) = 2cos(7r/(p + 1)), if we parametrize as usual 
c = 1 — Q/{p{p + 1)), and h = hi^2- P will play the role of the weight of loops in the 
study of loop models in the next sections. 

A small-g expansion of the amplitudes Aa can also be computed using the explicit 
expression of the associated boundary states (see Appendix B). To first order we have 



-4/11 



-c/48 



1 + 



A^{t) = 42'i2-4hi^-c/48+/i2/2 



(3-7p)2(-2+p) 



2p(l +p)(3 + p) 

9 3 

2 p 1 + p 



1 + 



q + ... 

45 



+ 



80 



l + 3p 



+ 



(97) 
(98) 



We checked agreement with the expansion of eq. (95) up to order q^^. 

Finally, to interpret the partition function above as that of the Potts model, we fix 
the normalization following [11], demanding that in the limit C — t- of an infinitely long 
rectangle, the normalized conformal blocks^ behave as 



AA«g"(l-C) = l + 0(C). 



(99) 



This holds because for a very long rectangle, one has in both cases a = 0, a = 2, the 
partition function of an horizontal strip with free boundary conditions on the sides. Both 



Redefining <p1 — > A" 0° and |0) — >■ a|0) we can fix opportunely say A^^j 
OPE structure constants give the behavior of eq. (99). 



"'^''^ and AJ'//'^ so that the 



-1/2 
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conformal blocks should give one in that limit. We have: 



1 



(100) 



1 
F20 




(101) 



4 The case of free theories 



We focus now on the description of rectangular amplitudes in free theories. This section 
completes the discussion of free theories presented in [3] with examples with different 
boundary conditions along the rectangle. 

4.1 The Laplacian of the rectangle 

We start our discussion of free theories by computing the partition function of (unrooted) 
spanning trees on a rectangle. This problem is described by a free fermionic field theory 
via the matrix-tree theorem, stating that the number of spanning trees on a graph G is 
given by the determinant of a minor of the Laplacian matrix A of G (see e.g. the review 



the zero eigenvalue. Then the number of spanning trees is given by a two point function 
in the fermionic theory. Neumann boundary conditions (A^) for the field translate to free 
for the spanning trees, while Dirichlet ones {D) impose that boundary sites are linked 
to a single external vertex ("wired" boundary conditions). 

We would like now to compute the Laplacian with either or D on the sides of a 
rectangle and take the continuum limit. This can be done by taking the tensor product of 
Laplacians on a segment, and from that one can extract the expressions for the universal 
part in the limit of infinite rectangle with fixed r = iL' /L. Some of these results were 



[35]): 




(102) 



where Ay is the matrix obtained from A by deleting the ith row and column, to eliminate 
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first given in [16] and one finds: 



(103) 



D N 



Vl 




(104) 



J_[D=iTT 

^ TV 



oc 



v{r/2) 
7?(r) 



(105) 




N 



oc 



r7(r/2)r/(2T) ' 




(106) 



In these results we did not keep track of proportionality coefficients, which are related 
to the OPE structure constants, and in principle different in each case^. The partition 
functions above are enough to determine all the other obtained by permuting boundary 
conditions. This is possible using the symmetry of describing the same partition function 
choosing imaginary time flowing either bottom-top or left-right (crossing symmetry of 
correlators of boundary fields), and modular covariance of these partition functions. For 
example from modular inversion one has: 



where we have used identities listed in Appendix A. Note in particular that the partition 
functions (103) and (106) are modular invariant. The most important physical feature 
exhibited by (103)-(106) is the dimension-full factor , with w the modular weight given 
by eq. (32). Indeed in the continuum the fermionic theory describing spanning trees is 
the logarithmic CFT of symplectic fermions with c = — 2, where the twist field changing 
boundary conditions N to D and vice versa, has conformal dimension hi^2 = — 1/8 [25]. 

Having discussed the case of c = —2 one can immediately find the result for c = 1, 
the free boson, by taking the power —1/2 of the results above since the twist field in this 
case has dimension 1/16 = (—1/2) (—1/8). The case of a free boson on a rectangle was 
discussed in detail in [22], and our results are in agreement with those findings. Further, 
in that work the authors discussed also the boundary states associated to different choices 
of boundary conditions. The results resemble those for the cylinder boundary states 
except that, as expected from the general discussion of section 2, for the rectangle only 

^ Moreover the free energies computed from the lattice Laplacian contain also an additional geometry- 
independent constant, which we cannot predict using CFT. For example in the case of D boundary 
conditions around the rectangle it is [16] Fo = ilog(4v^)- For a recent discussion of this issue in the 
case of the Ising model with free boundary conditions, see [36]. 




(107) 
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one set of oscillator modes an of the bosonic field is left. The boundary state describing 
the bottom of the rectangle is [22] 



= exp ( -eiet ^al] I^^^^J > (108) 

\ n>0 / 

where e^r,b equal to +1 (—1) means N (D) boundary conditions at left, right and bottom 
side respectively, and is the zero mode, which is the simple Fock vacuum except 

when ei = €r = +1, in which case is \p = 0) for ef, = 1 and \x = xq) for eb = — 1- 

4.2 The Ising model 

In [3, 23] the case of the Ising model with D boundary conditions on all sides of the 
rectangle was considered, corresponding to fixing Ising spins to the same state on the 
boundary. The result for the corresponding rectangle boundary state has the following 
coherent state form 

\B^) = exp j G„,„V'-m-i/2V'-n-i/2 I \0) , (109) 

\o<m<n I 

with \0) the vacuum annihilated by the positive modes of the Neveu- Schwartz fermions, 
V'p+i/2|0) = 0,p G Z>o, and Gm,n defined by 



^^^^'^^^ = ^(i^ 1-^ ^'''^ 

^ ^' \ z\Z2 j ^ ^ m,n=0 ^ ^ 

Here we will generalize this result to the case of different boundary conditions on the 
sides of the rectangle. This will be done by mapping correlators from the region with 
corners to the upper half plane, through the conformal map /(z) = (2 + 2;^^)/2, similar 
to the situation of figure 4. 

Let us consider first the upper half plane H with boundary conditions for Ising spins 
free / on {x < 1} U {x > 1} and fixed (say +) in between. This corresponds to inserting 
a BCC operator a^'^ at x = ±1. From the fusion rules of these BCC operators 

CT-^'+0f7+-^ = /-^'-^ + ^^'-^, (111) 

we know that the corresponding rectangle boundary state l-B^^) is composed by two 
basis states, one in the sector of the identity \ J and one in that of the energy \ 
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The numerical value of C is 



C 



C 



f,+J 

<J,a,ip 



C: 



fJJ 

o-.cr,I 



f,+,f 

cr.cr.I 



fJJ 



(113) 



as follows from using F-matrices for the Ising model (see e. g. [28]) after identification of 
boundary conditions with Virasoro representations f = a, + = I and normalizing two 
point functions to one. 

In order to evaluate the first contribution to the boundary state \B^j) we 

consider the two point correlation function of fermions in the upper half plane with 
these boundary conditions (see e. g. [1], eq. (22)) 



(V'(U'1)V'(W2 ))/+/,] 



1 



1 +Wl 



2{wi - W2) \-\l + W2 



1 + W2 

l + Wl 



1/2 



1/2 



1 -W2 
1 — Wl 

1 — Wl 
I -W2 



1/2 



+ 



1/2 



(114) 



where we obtained this result from the chiral correlator {'tlj{wi)'ip{w2)cr{l)(T{—l)) and 
omitting irrelevant factors. This result can be related to the boundary state l^jj) 
describing the geometry T> with / + / boundary conditions through 



{^PizMz2))f+f,v = {OmzMz2)\B+f) 



dwi dw2 
dzi dz2 



1/2 



{ll^{wi)tp{w2)) f 



(115) 



Following the discussion in [3], we look for a fermionic coherent state form of the bound- 
ary state 



I ' ) 

\tJ (J I 



4 '^^'^ : exp 



dz 



dz' 



i,{z)G{z,z')^{z')] : |0), 



(116) 



2iTT J 2i7r 

where the factor ^r'^'^" is introduced in order to obtain consistency with the definition 
of basis states (17). Substituting this in equation (115) one can easily determine the 
function G{z,z'). It is convenient to split expression (114) into the sum of two terms. 
Going through the calculations leads to the introduction of two functions 



and 



1 



4(^2 - Zi) 



A{Z2 - Zi) 



Z2 ' 



1 

^1 



7^ 

^2 



1 

Z1Z2 



(117) 



(118) 
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for which we define the expansions 



1 yr, ^-^ yl 



2ziZ2 



(119) 



-^1 •^Q 

m,n=0 ^ ^ 

SO that the boundary state gathers the contribution 



/.) = 4-2''- exp [^-«] V'-m-l/2V'-n-l/2 |0) 



(120) 



. 0<m<n 

where we defined G = + G^^), The first values of Gmn read 

Gqi = 

7 

Cos 



3 

'2 
3 



07 



Go5--Y^,Gi4---,G23--- 

75 3 55 63 

'128 ' " " 128 ' " "128 ' ^^'^ " "128 ' 



This leads to the first few terms in the expansion of the boundary state 



3 7 3 

1 + -V'-l/2V'-3/2 + ^V'-l/2V'-7/2 + ^^-3/2'0-5/2 + 



\0) 



:i2i) 



:i22) 



We now determine the other contribution \J'a) to \B^f) in addition to \J ^j). For doing 
that we use now the chiral correlator {'il){wo)ilj{wi)ilj{w2)(j{l)a{—l)) in [1] (eq. (66) of 
that reference) and send wq to infinity to obtain 



('0(oo)V'(u;i)V'(tf2))/+/,i 



wl+W^— W1W2 — 1 



{wi - W2){wi - ly/^iwj - ly/^ ■ 

Setting u = and u = we find the corresponding correlator in the z plane 

{v — u){l — uv) 



(123) 



Z2 ■ 



{'ip{oo)i;{zi)^{z2))f+f. 



V 



(l_„2)l/2(l_ ^2)1/2 
(1- 7x2)1/2(1 _ ^2)1/2 



+ 



{v — u){l — uv) 
We represent this expression alternatively as 



9 1 — ti^ 9 1 — 
V z + u 



1-w^ 



1-^2 



(124) 



(V'(oo)V'(zi)V(^2) : exp ( ^ ^ I ^^(z)G{z,z')i;{z') ) : ^(0)|O) 



2i7r / 2i7r 



(125) 
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This gives us the lengthy expression 

1 {V-U){1-UV) 1 (l-n2)V2(l_ ^2)1/2 

(^{Zi, Z2) - 2 (i_ ^2)1/2(1 _„2)l/2 + 2 (^;-u)(l-H 



' 2 1 - , 2 1 - 

^ 2 +^ 2 

1 — 1 — li^ 



where we used that 

(V^(00)V'(Z1)V'(Z2)V(0)) 

We finally expand as before 



+ ti^ — uv 



V — u 



G(zi,Z2) = Y 5] G„nn"^t;", 



(126) 
(127) 

(128) 



m,n>0 



from which we get the final result 



Icr (7/ 



^h^-2h^ exp 



G'mnV'-m+l/2V'-n+l/2 



0<m<ri 

In the end one gets the following expression 



m\o). 



(129) 



Icr cr/ 



1 



1 - ;7V'l/2V'-5/2 - g^l/2V'-9/2 - ^1/2 V"- 13/2- 



-V'-3/2V'-5/2 - ^W-3/2i^-9/2 + 



^-1/2|0) 



(130) 



for the second (odd) contribution to the boundary state. 

Taking overlap of these boundary states we can compute partition functions with 
different boundary conditions of Ising spins, and their agreement with four point func- 
tions of BCC operators along the general scheme of section 3, can be readily checked. 
We will not present here the details of this check, but support our results with lattice 
computations for the Ising model instead. 

4.2.1 Lattice Ising model 

In this section we will compute the scalar products of the discretization of the boundary 
state \Bj-f) with excited states in the Ising Hamiltonian. We consider the Hamiltonian 
limit of the 2D critical Ising model on a strip with free/free boundary conditions: 



H 



1 



L-l 



\j=i i=i 



(131) 
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Now recall that the Hamiltonian (131) is related to the 2D Ising model in the basis. 
Then to make the contact with the 2D model, we should rotate 90 degrees clockwise 
the spins of the chain in the x — z plane, that is, we define the "fixed" boundary state 
\Bjf)L as 



1 

{Mf=t,4} 

We diagonalize the chain [29] to obtain the eigenvectors \k)L, and compute numeri- 
cally the scalar products L{B~j'j\k) l. This task is simplified by using the free fermionic 
nature of the problem, in a way similar to what we did for free boundary conditions in 
[3]. This time, however, the computations are more complex and we have restricted the 
analysis to small system sizes (up to L = 30 sites). Before presenting the results, we 
recall that due to the anomaly at the corners (see eq. (32)), scalar products between a 
lattice rectangular boundary state \B)l and the (normalized) fc-th excited state \k)L of 
the strip Hamiltonian, are expected to scale as [3]: 

-\og{L{B\k)L) = aoL + ai\ogL + a2 + j- + ^ + o{^j^ , (133) 

with 

ai = 2{hi + hr)-^. (134) 
o 

The value of the CFT scalar product {B\k) is then extracted from 

a2 = 7-log{{B\k)), (135) 

7 determined normalizing {B\0) = 1. See section 5.2 below for a systematic study of 
these coefficients. For computing {B\k), we actually fit 

, / L{B\k)L\ as 04 as 

We present the data obtained in table 1. The value ai = 2 x 2/16 — 1/16 = 3/16 is 
found with good accuracy only for the first lower excited states. The continuum limit 
of the excited states of the chain is associated to a CFT state obtained by the action 
of negative fermionic modes on the vacuum. It is important to distinguish the two 
sectors, even and odd number of fermions, corresponding respectively to descendants of 
the identity I and of the energy -0 in the CFT. We replace 7 in (135) by 7 + 7 if is 
a descendant oi ip. 7 is written as minus the logarithm of a prefactor C multiplying the 
states in the sector of ip. Looking at {B^j-\1) (where |1) = V'-i/2|0)) determines 

C = 1.41422 ±0.00002. (137) 
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This agrees very well with the prediction from the general theory eq. (113) (the factor 
4-/1^ simplifies with 4^^!^ -2/1(7^4-2/1,7 f^om the normalization of the basis states la^a) and 
lo-^o-))) providing a lattice measurement of boundary OPE coefficients. In table 1 we 
list the scalar products obtained. The agreement with the CFT prediction is good, 
validating our analytical computation. 







numerics 


CFT 




0) 





1 


1 




1) 


1/2 


1 


1 


Hf 


2) 


3/2 








Hf 


3) 


2 


1.499999 ± 0.000020 


3/2 = 1.5 


Hf 


4) 


5/2 


0.500079 ± 0.000012 


1/2 = 0.5 


Hf 


5) 


3 








Hf 


6) 


7/2 








Hf 


7) 


4 


0.875178 ± 0.000020 


7/8 = 0.875 


Hf 


8) 


4 


0.374991 ± 0.000005 


3/8 = 0.375 




9) 


9/2 


0.375315 ±0.000031 


3/8 = 0.375 


(S+llO) 


9/2 


1.125002 ±0.000020 


9/8 = 1.125 



Table 1: Scalar product of l-B^^) and \k) from finite size scaling (numerics) and compar- 
ison with CFT prediction, h/. is the conformal dimension of the field \k). Chains of size 
up to L = 30 sites are used. 

5 Loop models 

5.1 BCFT of loop models 

The lattice model we consider now is the dense loop model based on the Temperley- 
Lieb algebra underlying the Potts model [30] . We consider a system of L (even or odd) 
strands, with free (reflecting) boundary conditions at both boundaries. At the critical 
point, the anisotropic version of the model is defined by the Hamiltonian: 

L-2 

H = -Y,ei, (138) 

1=0 

where Cj are the TL generators. This Hamiltonian acts on link states, patterns of con- 
nectivities of sites, and has a triangular block form, with j = 0,2, ... ,L (L even) or 
j = 1,3, . . . , L (L odd), the number of through lines, indexing each block [30]. 

The continuum limit of this loop model, when we parametrize the loop weight as 
/3 = 2cos(7r/(p ±1)) (p > 1 is a real parameter which we keep generic here), is a CFT 
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with central charge c = 1 — ()/{p{p + 1)). The primary boundary fields of the resulting 
(quasi-rational) boundary CFT are (pj, j € N (see e. g. [33]). is interpreted as the 
operator creating j lines in the continuum limit of the loop model, and corresponds to the 
irreducible Virasoro representation with highest weight <Ai,i+j (Kac table notation), with 
conformal dimension = '^^^^^jy- The primaries <f>j generically fuse as spin-(j/2) 

su(2) representations: 

VixVj = Y,Vk, I = {\i-j\,\i-j\+2,...,i+j}. (139) 

The allowed boundary conditions one can put on the strip are given by the labels j 
of Virasoro representations. Setting i and j on the boundaries will then allow the 
propagation of only the sectors given by the fusion of i and j [8]. It is possible to 
manufacture microscopic boundary conditions for the lattice model to reproduce this 
behavior [32] (and see below). Note that free (reflecting) boundary conditions on the 
lattice, without restrictions on the number of through lines, correspond to setting on 
both boundaries the label j — )• oo. The BCC operator between boundaries j and j' is 
given by 4''lj_j/\ i the operator of smallest conformal weight in the spectrum of the transfer 
matrix for a strip with j and j' on the sides. Fusion rules of these boundary operators 
are only a subset of those for free boundaries, eq. (139) (because some boundary OPE 
structure constants vanish): 

seX' 

where X' = ({\i - k\,\i - k\ + 2, . . . ,i + k} n {\\i - j\ - \j - k\\,\\i - j\ ~ \j - k\\ + 2, . . . ,\i ~ j\ + \j - k\}). This 
restriction comes from the compatibility of the number of lines inserted by the fields in 
the corners and that which is allowed by the boundary conditions. We depict in figure 
6 the geometric interpretation of the two-point function of BCC operators as the con- 
traction of \i — j\ lines originating from the point in which the boundary condition is 
changed from i to j and terminating at the point in which we change from j to i. More 
details about this geometric interpretation of the fields will be given in the next section, 
where we also present numerical checks. 
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Xi X2 



Figure 6: Geometric interpretation of the two-point function of BCC operators. 
5.2 Numerics 

We now want to understand which states on the lattice are described by our basis 
boundary states introduced in section 2. Define the normalized link state \i^j)L as: 



U J/ 




with 



Here, 



ni + ric = i ] rir + ric = j ; ni + rir 



(141) 



(142) 



lC s\i^j)L is the lattice norm of the state. It is obtained by using the loop scalar 



product of link patterns [14], which turns out to converge nicely to the Virasoro bilinear 
form in the continuum [14]. We call the s = ni + nr lines not paired, through lines. The 
link state we consider is projected onto the sector with exactly s through lines, and each 
time the through lines are contracted we get zero. The ric lines inserted at the left corner 
are paired with the ric ones inserted at the right corner, and during the evolution, these 
lines can propagate through the system or be contracted. Now we claim that this state 
flows in the continuum limit to the following basis state: 



(FP) 

>L — > a 



(143) 



The limit to be taken is L — )• oo with z, j, s fixed, and we keep only the finite part (noted as 
(FP)) which is supposed to be described by the CFT, as will be discussed more precisely 
below. This identification corroborates the interpretation of the boundary operators 
and as inserting respectively j = n,. + ric lines at the right corner and i = ni+ric lines 
at the left corner. We remark again that the paired ric lines on the left corner can be 
contracted with the ric lines on the right corner during the evolution, but the remaining 
s = rij. + ni are forced to propagate in the system. Our discretization of the basis states 
is a natural generalization of that for the case of no line insertion discussed in [3] . 
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We support our claim with the numerical results (labeled num) of table 2, obtained 
for fugacity of loops /? = -v/2 and to be compared with the CFT results (labeled theo), 
corresponding to central charge c = 1/2. Denote by \k)L the k-th eigenvector (normalized 
according to the loop scalar product) of the Hamiltonian restricted to the sector with s 
through lines. As discussed above, the scalar product of a rectangular boundary state 
and has the form (133), and now 7 = — log (a), a := a4'^''~^'~^^ , see eq. (17). In 
the data presented below a is determined by giving the expected value of ai (ai theo) 
as input in the fit. 

The scaling limit \k) of \k)L used to determine {B\k) (theo), are for generic k given 
by a combination of Virasoro descendants at a level determined by the energy of The 
state |0) is clearly identified with the primary \(j)j) for each sector. Then recall that the 
continuum theory has field content given generically by the quotient of Virasoro Verma 
modules T/i^i+s/Vi^„i_s (with the usual Kac table notation for Verma modules) [33]. 
In the sector s = 0, we identify then (due to normalization) |1) = -y/2/cL_2|0), |2) = 
yT72cL_3|0). When s = 1, we have |1) = l/VW^-il'/'i), |2) = l/v^4/ii +c/2L_2|0i). 
For s = 2, |1) = l/y/2h^L_i\(p2), and |2) = 2/3L_2|02) at c = 1/2 (note that this last is 
a priori an unknown combination of Virasoro descendants at level two, which for c = 1/2 
is fixed using that |(/>2) is also degenerate at level two, since /ii^3 = /i2,i)- Scalar products 
{B\k) (theo) in table 2 are then computed using the above formulas for |A;)'s and the 
explicit form of |-B)'s, eq. (17) (see Appendix B for and eq.(86) for Iq^^^) and 

1 2^0))' specialized at p = 3. 





L 


ai num 


a\ theo 




(B\l) num 


(fl|l) theo 


(B|2) num 


(S|2) theo 


ll 1. 


)l 


0.18712(307) 


0.1875 


1.27053(16) 


1.50004(21) 


1.5 








1 1 ■ 
lo 1 . 


)l 


0.06278(131) 


0.0625 


1.08322(6) 


0.70704(5) 


~ 0.70711 


0.35360(3) 


~ 0.35355 


1 2 ■ 

I2 0. 


)l 


0.93778(865) 


0.9375 


4.33940(152) 


2.00057(74) 


-2 


2.50616(105) 


2.5 


I 2 ■ 

II 1. 


)l 


0.18938(76) 


0.1875 


1.79691(8) 








0.50115(15) 


0.5 



Table 2: Lattice scalar products for loop fugacity /? = ^/2 and comparison with the 
results for the c = 1/2 CFT. Extrapolations are determined by fitting values for L = 
10 24. d|^. = af^A^'^-^'^-^K 

Finally, we note that the state defined in (141) is of course not the only microscopic 
state giving rise to 1^^^) in the continuum limit. Any similar state obtained via local 
modifications (e. g. by allowing a finite number of arches between the through lines, 
etc.) would obey the same property, albeit with a different value of the coefficient af^. 
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5.3 Geometric interpretation of conformal blocks 

Up to now we have given numerical support to our claim that the rectangle basis states 
I j defined in eq. (17) can be interpreted as the continuum limit of microscopic boundary 
states with a given pattern of connectivities. Here we want to push this discussion further 
and associate a geometric picture in terms of loop configurations to conformal blocks 
building the correlator of BCC operators. According to the previous discussion, we have 
the identification: 

LCi|T^'|/fc)L a^atiA'^^^'^'^ir) . (144) 

We depict in figure 7 two configurations contributing to the lattice amplitude l (^2^1^^' 12^2) l 
and one to the amplitude l{^4^\T^' {2^2) l, drawing only lines inserted at the corners. 

i = 2- 

s = 2: 

i = 2- 




Figure 7: Configurations contributing to the lattice amplitudes l{'^s^\'^^' {2^2) l for s = 
2,4. 

We see that AI''^''''\t) is then associated to the continuum limit of lattice amplitudes 
where we insert i lines in top left corner, / in top right, j in the bottom left and k in 
bottom right, s of the lines (the most exterior ones) are forced to propagate. The other 
lines can — but are not obliged to — be contracted. Only if s = i + I = j + k (s = 4 in 
the example of figure 7), the amplitude is associated to configurations where all lines 
inserted at the corners are forced to propagate. For example, the CFT partition function 
of figure 2 will be Z (0 • j+j) oc L'^Al'^j''^. The only term in the amplitude sensible to the 
changing of the number of through lines (associated to different fusion sectors) is the 
conformal block. Then we are led to the geometrical identification of conformal blocks 
which we schematically represent as 
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where the diagram on the right has a fixed number s of through hues and the bubble 
stays for the possibihty of contracting or not the other hues inserted at the corner. 

We recah that in general the partition function Z is given by a sum of amplitudes from 
the relation (31), and its geometric interpretation follows from that of the amplitude. 
We will comment more on that in section 5.5, where we discuss the universal character 
of the coefficients af^, and the precise relations between lattice amplitudes and CFT 
partition functions (30). 

5.3.1 Probabilistic interpretation 

We can easily obtain interesting results from the discussion above as follows. We spe- 
cialize to the case of one line insertion at every corner, i = j = k = l = l. In this case, 
the amplitudes are formally equivalent to that for the Potts model with insertions of 

f F 

4)12 changing from free / boundary conditions to fixed F, and were computed explicitly 
in section 3.2.2. The partition function 



is unambiguously the (universal part of the) continuum limit of a loop model where 
the two lines inserted at left and right bottom corners constrained to propagate in the 
vertical direction. Consequently 



is then the partition function of a loop model with the two lines inserted at left bottom 
and left top corners constrained to propagate in the horizontal direction. Then one can 
interpret the conformal blocks Q"" as in the top of figure 8. The geometric interpretation 
of every other amplitude involving boundary states with insertions at the corner will 
follow from the decomposition in terms of ^2(7") and A2{—^/t). For example for Aq{t), 
we can use the hyper geometric identity ^'^(l — C) = 1/-^20^^(C) ~ -f22/-F20^^(l — C) oc 
g^iC) + l//3a^(l - C), see the bottom of figure 8. 

We have seen that the two elementary events which can happen when inserting a line 
in each corner are described by Q'^{1 — C) and ^^(C)- Then, the probability of having 
these two lines fiowing vertically is: 



Z (0 1 2) (r) = M^r'^\T)ef^{T)G\l - C) oc L"'^2(r) 



Z{ili) {r)=M^L-'r^--/\T)ef\T)g\Q{T)) cc (L')"'^2(-1/t) 






+ 




Z(o{2)(t) + Z(iOi)(t) gi(l-C) + gi(C) (145) 
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Figure 8: Geometric interpretation of the conformal block and G^. Thick hnes repre- 
sent the non-contractible hnes created by the one-leg operators in the corners propagating 
or not through the system (time is flowing upwards). 



The probability that the lines propagate horizontally is clearly 1 — P{t). Note that 
whatever normalization of the partition function drops out since it is common to both 
partition functions. 

We first consider the case of percolation, for which /? = 1, corresponding to c = h = 0. 
We have: 

p (^\_ (i-C)^2i^i(i|;|;i-C) , . . 

(1-C)'.F,(1 |;|;l-<)+d.F,(l,|;J;C) '"'^ 
- „_04,F,fl,|;il-c) (147) 



r(|)r(i)' ^' V3'3'3 
rQ)r(l) 



(1 - C)^ (i - ^(1 - C) + |(i - 0' + o ((1 - C)')) . 

(148) 



This quantity is precisely the crossing probability computed by Cardy in [9, 11], that 
is, the probability that there is at least one Fortuin-Kasteleyn (FK) crossing cluster 
spanning the rectangle in the vertical direction. We now comment on this relation. 
Microscopically, the configurations of loops contributing to the numerator of our formula 
are those where the two lines inserted at the corners propagate from bottom to top. In 
terms of FK clusters (which are encircled by the loops), this corresponds to counting 
cluster configurations with "wired" boundary conditions on top and bottom rows, and 
to the constraint that there is at least one spanning cluster crossing the rectangle, see 
figure 9. Call such ensemble of cluster configurations Syj. 

Instead the configurations counted by Cardy's formula, which we denote by S, are 
simply those with at least one spanning cluster, without imposing wired boundary con- 
ditions. We now show that 



#(5) = 2^ X 2^ X #(5^) , 



(149) 
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T = L' 



r = 




Figure 9: A configuration of loops (black) and the corresponding configuration of 
Fortuin-Kasteleyn clusters (gray). The through lines inserted at the bottom corners 
are forced to propagate to the top, and wired boundary conditions (all bonds open) 
are imposed on the links of the bottom and top rows of the lattice where the Fortuin- 
Kasteleyn clusters live, r = and t = L' are initial and final discrete imaginary times 
of the transfer matrix evolution in the vertical direction. 

where #(-E) is the number of elements of the set E. First note that changing the status 
(open or closed) of links in the bottom and top rows of the lattice where FK clusters live, 
transforms a configuration s £ S into another s' G S. Then we can group the spanning 
configurations in classes differing by only the status of bottom and top links, and for 
each class we choose the representative with all the bottom and top links open. Since 
there are 2^ x 2^ (first factor for the possible status of top links and second for those 
of bottom links) elements in each class, eq. (149) follows. Now note that also the set 
of all possible (not only spanning) cluster configurations can be divided in classes with 
elements differing by the status of top and bottom links as before. Then, if we compute 
the crossing probability by dividing the number of spanning configurations #(5) by the 
number of all possible configurations, the multiplicative factor 2^ x 2^ drops out, and the 
equality of our formula with Cardy's one is established. The above argument does not 
work anymore when we have generic weight of clusters /?, since spanning configurations 
cannot be regrouped then due to the different powers of (3 weighting each configuration. 
Crossing probabilities for generic fugacity /3 have been computed using SLE methods, 
see [2]. 

5.4 Logarithmic cases 

We have anticipated in section 3.2.2 that for certain (logarithmic) values of the central 
charge the conformal blocks for insertions at the corners, (93)-(94), are ill-defined. 
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Indeed, more fundamentally fusion of the fields at these points cannot be decomposed 
onto irreducible Virasoro modules [20,21]. For the operators degenerate at level two we 
are considering, this happens for c = —2 and (pi = (j)i^2 h = —1/8 (dense polymers), 
and for c = and (pi = </>2,i, h = 5/8 (dilute polymers). In these cases, we have the 
following behavior of conformal blocks (e = p — 1 in (150), e = p — 2 in (151), where we 
parametrized c = 1 — 6/(p(p + 1)) 

limgO(c) = -K{C) ; i^gHC) = -K{C) ; (150) 

e-)-o vr e-i>o vr 

1 

Q'iC) = TTrS2{C) + 0{e) ; lim ^^(c) = 5^(0 • (151) 

Ibe e^o 

is the complete elliptic integral of the first kind of parameter C, and 

.2^/13 



S2{0--=C2Fi[--,--?,-Cj . (152) 

For dense polymers the degeneracy Q^{C,) ~^ ^^(C) arises because the roots of the indicial 
equation (92) coincide, so that one has to replace one conformal block with the missing 
solution of the differential equation which generally is given by a derivation procedure 
w.r.t. the root, and here can be written as (e = p — 1) 

,__^g^(0-g°(C)^ logW;,-(C). (153) 

e->0 e vr 

In terms of fusion rules this is stated as the degeneracy of the two fields </>o and (j)2 , which 
are replaced by the identity and its logarithmic partner under the Virasoro algebra [20] . 
For dilute polymers, here has to be replaced by the logarithmic conformal block, 
which can be written as {e = p — 2) 

" " ^ " ("^' ^' ^' ^ " ^^^^^ 

= ^2(1 - C) , (155) 

where Fij are the fusing matrices of four (pi fields at generic central charge, eq. (96). 
Again this is stated as the degeneracy of the stress tensor and the field (p2, which are 
organized in an indecomposable module under the Virasoro algebra [21]. 

The presence of a divergence is clearly a difficulty for the interpretation of the parti- 
tion function as that of a statistical model. We consider the partition functions Z (o ]^ a) , 
a = 0, 2 computed in section 3.2.2, for which we have the geometrical interpretation 
discussed in the previous section. In particular recall that Z (o } 2) is given by config- 
urations of loops where we constrain the two lines inserted at the bottom corners to 



5 LOOP MODELS 



35 



propagate without being contracted. We assume now that the geometric interpretation 
of boundary states goes through also in the dilute case, where the value of of h in the 
formulas is /i2,i. We expect physicahy that these partition functions should be well de- 
fined for every value of loop fugacity /?, also at /3 = 0, which is the case we are after 
for polymers. The solution to this is simply hidden in the normalization of the partition 
function, eq. (95). In the case of dense polymers = N"^ = 0. This does not come as 
a surprise, and should be traced back to the presence of a single loop on the rectangle. 
Factorizing this trivial zero, one then gets that for dense polymers: 

Z (o i 2) (/i = -1/8) = Z (0 } 0) (/i = -1/8) = ^/Lrie^^-K{l - Q (156) 

= , (157) 

where we used the identity 2K{1 — C,{'t)) = 7r^3(T)^. As a highly not-trivial check of the 
above result we can compare with the expression coming from the determinant of the 
Laplacian with different boundary conditions on opposite edges, see eq. (104). 
For dilute polymers one has 

AA0^|e + O(e2) (158) 

AA^^^7r + 0(e). (159) 

Note that the indecomposability parameter h = 5/6 [21], characterizing the Jordan 
cell of the stress tensor for dilute polymers, is present explicitly in the coefficients, 
h^/b = 15/32. Finally we predict again degeneracy of both partition functions 

Z{o\2){h = 5/8) = Z{o\o){h = 5/8) = L-^9^^^S2{1 - C) ■ (160) 

The degeneracy we have observed could be understood more fundamentally as the gluing 
of the two sectors 0o and ^2 in an indecomposable module corresponding to a single 
boundary condition. 

It is of interest now to see what is the behavior of probabilities (145) at the loga- 
rithmic points discussed above, c = —2 dense and c = dilute. The outcome is that 
probabilities are well defined but their expansion shows logarithmic terms. For expand- 
ing around ( = I, the case of a tall rectangle, we use 

GHO = F2oQ\l - C) + F22Q\l - . (161) 

As /i — )• — 1/8 (dense case) we see from formulas (96) that F20 — )• 00 and F22 — )• —00, but 
-^20 + -^22 — log(16)/7r, and recall that the conformal blocks are degenerate, eq. (150). 
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The probability reads: 

Pc=-2,dense(C) = ^(^^^J) + ]^(^) (162) 

+ , "^'"5^ ..2 +0(a-C)^)- (163) 



- - log (y ) 2 - log (y 

When instead /i — )• 5/8 (dilute case) we have F20 — )• b/h'^ = 32/15 and F22 — )• —00. 
In this case the divergence of Q^, eq. (151), compensates that of the OPE structure 
constant, and the result is finite: 

Pc=o,dnute(C) = sJ-0~+%{0 ^'^^^ 
= ^(1-0^ + ^(1-0^ (165) 

757r f 12 log (^) + 127r - l) 
+ ^ ^ ^(1 -Cr + O ((1 - 0') • (166) 

We are not aware of a previous derivation of the above probability formulas. Note 
that since P(l — Q = 1 — P{C), the same series plus the constant term holds for a very 
long rectangle (^ — > 0) replacing 1 — by C- We have already seen how rectangular am- 
plitudes contain the OPE structure constants, and now we see that for logarithmic CFTs 
the indecomposability parameters (contained in the i^-matrices) show up in geometric 
observables. 

Note also that expressing the anharmonic ratio in terms of the modular parameter 
q through eq. (44) 

log(l - C) = log(16) + log(Vg) -8^+... (167) 

integer powers of r enter explicitly the expression of the probability. This feature is a 
consequence of the Jordan form of the transfer matrix in these cases. 

More generally one could replace the ill-behaved conformal blocks in the logarith- 
mic cases by a combination of them which has probabilistic interpretation and which is 
expected to be well defined even if the conformal blocks themselves are singular or de- 
generate. Unfortunately writing down explicitly such a probabilistic basis in the space of 
conformal blocks does not seem to be simple for higher number of line insertions. Indeed 
we do not know how to associate precisely a conformal block to a given propagation of 
lines (see discussion at the beginning of section 5.3), but more importantly, we do not 
know how to fix the normalization. 



5 LOOP MODELS 



37 



5.5 OPE coefficients and lattice crossing symmetry 



Given a CFT rectangle boundary state = s ^fjs \J ^ssS I i ) ' discretiza- 

tion of basis states introduced above in eq. (141), we can define the more general lattice 
boundary state as combination of lattice basis states 

\BL)L = Y.D^Js\^'J)L■ (168) 

s 

We expect then that \B^^)l ^—l if 

1 



ij 

We now question the universal character of a. As the crossing symmetry of correlation 
functions (or rectangular amplitudes) is a powerful tool for constraining OPE structure 
constants in the continuum theory, we will see that also lattice crossing symmetry allows 
to derive constraints for the coefficients L>'s appearing in eq. (168) and to infer about 
the universality of the coefficients a's. 

Let us see concretely how to implement lattice crossing symmetry for the simplest 
non-trivial case of i = j = 1, corresponding to boundary fields inserting one line in each 
corner. For this situation all the possible CFT rectangle boundary states one has to 
consider are: 



-^i,i+2/ = ^1,1,2 V'^220 111)' (I'LI) 



I R2±l\ _ /^iii i \ , ^j,«±l,j /^iii I 2 \ 

/ — <-^l,l,0 Y '-^OOoll 1/ + ^1,1,2 Y<-^220lll/' \^'^) 

where boundary labels G N. The discretization of the above basis boundary states follows 
from the discussion at the beginning of section 5.2 : 

^''^«?ilA), (172) 



(FP) „2 I 2 




pL/4-1/2 U • • • U 

Accordingly we introduce the discretizations 



«nlii)- (173) 



^i,i+l,i+2 

\BI^+2)l = ^L/4-1/2 I U • • • U I ' (^'''^) 



T-)i,i±l,i j-)i,i±l,i 

im±i\ fAio_ i_rli2_ (^7^] 

i^i.i )L- or. /A vu ••• u; +^L/4-i/2 I u ••• u I ■ ^^''^^ 
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L = 4 



Figure 10: A lattice of size (L = 4) x (L' = 6). At each node we have the transfer matrix 
depicted on the right, equal to identity plus TL generators. 

Now consider an L x L' lattice as in figure 10, where each vertex corresponds to a 
node transfer matrix equal to the identity plus the Temper ley Lieb generator. 

A partition function on this lattice is specified by giving an initial and a final state, 
and computing the loop scalar product of the final state with the sum over all possible 
choices of node operators acting on the initial state. The universal part of the partition 
function in the continuum limit is a partition function in boundary CFT. The conformal 
boundary conditions are imposed by the choice of boundary states according to the 
identification of their continuum limit as discussed in section 5.2. Clearly one can decide 
to consider time fiowing from bottom to top or from left to right; the equality of these two 
descriptions for each configuration is a constraint that we call lattice crossing symmetry 
and which allows for the determination of (ratios of) the coefficients -D's. See figure 11 
for a picture of the resulting equation. 



+ 
+1 



Figure 11: Equation imposing lattice crossing symmetry. 
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Consider first the case of boundary conditions i,i + — 1 around tfic rectangle. 
Taking tlie identity operator acting in tfie bottom-top direction (or TL generators in the 
left-right direction) at every vertex, gives the equation 

L'/2-l 

l{B^'\BI^')l = P'^^'-\'{BI^,,.^\ n e2^+i\Bl^,,.^L' = 0, (176) 

i=0 

leading to the first relation: 

Then consider the case with i,i + + 1 around the rectangle. Again taking identity 
in the bottom-top direction at every vertex gives 

L'/2-l 

L{Bi^\Bl^)L = P'''^-\'{Bl^i,+^ n e2.+i|i3i+i,+i)L', (178) 

i=Q 

yielding: 

' + (p'lit^'') ' = ^^/2-^'/2^ (i^no ''''^') ' • (179) 

Now, taking instead contractions in the bottom-top direction (identity for left-right 
direction) at every vertex gives 

L/2-1 

-Q e2,+i|Bi+i)L = L'(5:+i,,+iis:+M+i)L', (180) 

i=0 

^ ^ ^V/2-L/2^ [D'^t''")^ . (181) 

(182) 



yielding: 
Defining 



-^110 



these equations imply: 



Ci = J ^ ^-1, Ci = -^. (183) 

These are recursion relations for the coefficients C with initial condition C^J^ = 0, and are 
independent of the system sizes L, L' . The solution of the recursion can be expressed in 
terms of Chebyshev polynomials of the second kind Un{x) as: 



(184) 
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The first few coefficients read explicitly 



Ci = ./^^; Cl = -^^, (185) 



1 



Cl = VW^; Cl = --j==, (186) 
. + 1 o I /32 - 1 

We now claim that these expressions are not restricted to the two cases solved above and 
indeed solve every possible constraint one can write down from switching the left-right 
and bottom-top descriptions of lattice partition functions. Using (169) and setting i = I, 
we have that the ratio a^^/a^i is: 




^2 1 ^101 /TnTT 

"ll ^ ^ '-^112 V '-^220 /-,oo\ 

Note that crossing symmetry only allows to determine a's up to a common constant. 

As already remarked, if we had started with lattice boundary states differing from 
those we have used by local modifications (e.g. by allowing a finite number of arches 
between the through lines, etc.), we would expect that the continuum limit would be 
the same, but with different a's. For example, we could consider instead of eq. (172): 

(189) 



If we start with different boundary states, the equations resulting from lattice cross- 
ing symmetry would be modified in general (e.g. if we use the boundary state above, 
eq. (179) will have an extra factor of in the r. h. s. , similarly eq. (181)). The ratio 
of coefficients for boundary states with local modifications afj^/a^i should always be 
independent of system size and related to /3 and OPE constants, but in general different 
from (188). Since we expect more generally that afj changes under a local modification 
of the microscopic state (e.g. by extra factors of /3 in the example above), it cannot be 
predicted by the CFT. Note that the "physical" states \B^^)l defined in eq. (168) are 
normalized in order to obtain in the continuum the CFT states \B^^), and the CFT 
partition functions (30), when taking overlaps. 

The expression (188) can be simplified further using eq. (66) for i = j = k = l = l, 
a = c = 0, b = d = 1, r = 2: 



(Cno)''C'oooi^02 \ \ + (Cll2)''C'220-^22 \ l -0, (190) 



1 1 
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and the explicit results (96) yield 



a 



11 



/3r 



8hi 
3 



J r I 



8fci 
3 



(191) 



^ (/32 _ i)r (1 _ 4|i j r(-4/.i) 

In table 3 we check these results against numerical values obtained for p = 3, 4, 5, 
with p parametrizing the fugacity of loops as /? = 2 cos(7r/(p + 1)). Large sizes which we 
did not study would be needed to get more precise results, but nonetheless we find good 
agreement in these cases. 



p 


num 


theo 


3 


0.70715(9) 


~ 0.70711 


4 


0.68506(32) 


~ 0.69267 


5 


0.66671(66) 


~ 0.68736 



Table 3: Comparison of the ratio afi/aii determined numerical (num) by fitting values 
for L = 10 — )■ 24 and theoretically (theo) using eq. (191) for different values of p, 
parametrizing the weight of loops f3 = 2cos(7r/(p +1)). 

Now that we have determined the ratio a^^/a^x we can invert the reasoning and 
determine the OPE coefficients in terms of it from eq. (184): 



^112 



•-"220 



^ijiitlji 
*-^110 



•-"000 



r 



8fei 
3 



4fei 
3 



U, 



±1 \ 1 



(192) 



In particular the rectangle partition functions for 0i insertion in each corner can be 
explicitly written as (the functions ^'s are defined in section 3.2.2) 



i + l 



+ 2 = AL^/^-^'''7]-^/^el^^'T^{i - z) 



i + l 
i + l 



i + l 



BiL'/^-^^^rj-'/^el^''' -F°(l - z) +F( 



/3 



f/.-i ( I 



02 



(193) 



T\l-z) 



(194) 

The constant of proportionality Ai, Bi can be read off from the general formula (31) and 
are given in terms of boundary OPE coefficients. We note that for the CFT describing 
loop models we are considering here, these coefficients are known and can be obtained 
by analytic continuation of the F-matrices for minimal models A4.{p,p + 1) [19]. 
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6 Conclusion 

In this work we have constructed CFT boundary states encodmg rectangular geometries 
with different boundary conditions, and elucidated their relation with correlators of 
boundary conditions changing operators, completing the investigation initiated in [3]. 
We have developed in particular the geometrical interpretation in terms of lattice loop 
models, and we have derived new formulas for probabilities of self-avoinding walks. 
Remarkably, these formulas contain the indecomposability parameters of logarithmic 
CFTs in the expansion of the amplitudes. 

In conclusion, the formalism developed in [3] and the present paper gives access to 
large classes of loop models partition functions on the rectangle. A further generalization 
should allow one, for instance, to obtain closed formulas for the average conductance of 
rectangles at the higher plateau transitions in the spin quantum Hall effect, along the 
lines of [4, 5]. 
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A Useful formulas 

We collect in this appendix definitions and useful properties of special functions appear- 
ing in the main part of the manuscript. 

The theta functions we use are defined as: 



with q = e^'^*'^, r = L' /L. Denote by K the complete elliptic integral of the first kind 
with modulus k and K' := K{1 — k'^) its complementary, then r = iK' /{2K). One has 




(195) 



(196) 



(197) 



(198) 




K = 
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Modular properties of these theta functions are: 



rii-l/r) = ^/^riir) (200) 
92(-l/r) = V^94{t) (201) 



93(-l/r) = V^e^iij) (202) 



^4(-l/r) = ^/^e2(r), (203) 

r,(r + 1) = e*"/i%(r) (204) 

02(r + l) = e*"/402(T) (205) 

e^{T + l) = e^{T) (206) 

^4(r + l) = ^3(T). (207) 

B Boundary states for two one-leg insertions 

We give here cxphcit expressions of basis boundary state \^-^) and l^^^^) for generic 
c = 1 — 6/ {p{p + 1)) up to level 8: 



42"! \^°^) = |0) + (7 - ^) i_2|0> 
V P + 3/ 



/ 40 24 1\ / 200 48 19\ , 

+ i-4 + + — 

V3p + 5 p + 3 2/ V9p+16 p + 3 6/ 



j|0> 



128(p + l)p=^ 2 (p(p(709p + 769) - 1665) + 315) 
-Zv_3|u; H ; 



3{p + 3)(3p + 5)(5p+7) ' 6(p + 3)(3p+5)(5p + 7) 

(p(667 - p(103p + 331)) - 105) 3 64(2p - l)p 

6(p + 3)(3p + 5)(5p + 7) 3(3p + 5){5p + 7) 

(p-9)(p(p(709p+ 481) - 1761) +315) 64p (p (ep^ + p - le) + 9) 

i/_Q (J) + -L/_f!ij_2 U) 

12(p + 3)(3p + B)(5p + 7)(7p + 9) (p + 3)(3p + 5)(5p + 7)(7p + 9) 

80 196 324 24 1 \ , 

+ + + t: ^-4^1210) 



(208) 



3p + 5 5p + 7 7p + 9 p + 3 
800 1372 2916 4 11 



9(3p + 5) 3(5p + 7) 7(7p + 9) p + 3 504 
256p^(p+ 1) 



^--210) 



3(3p+5)(5p+7)(7p + 9) 



i-5i-3|0) 



14896 ... , -, , , 



9{3p + 5) 15(5p + 7) 7(7p + 9) p + 3 1260 
128(p-9)p=(p+l) 



3(p + 3)(3p + 5)(5p + 7)(7p + 9) 



l3L_2|0) + . 
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», », 9 5p - 1 2{p +1) „ 

3p + 1 3p + 1 

34p3 + 21p2 + 44p _ 3 16{p - 2)p(p+ 1) 

-^ — 4 02) — — 1 02/ 

60p3 + 86p2 + 40p + 6 (2p + l){3p + l){5p + 3) 

p=^+34p- 3 2 20p3 + 58p2 + 44p + 6 , ^ 2{p + l)^ (2p + 3) , 

30p2+28p + 6 30p3 + 43p2 + 20p + 3 ' 30p3 + 43p2 + 20p + 3 

64p(p(p(p(43p + 96) + 44) - 21) + 18) 



9(2p + l)(3p + l)(3p + 2)(5p + 3){7p + 5) 
16p(p + l)(p(p(17p - 9) - 48) + 252) 



-t-e 102) 
-t-5,-£'-ll02) 



9(2p + l)(3p + l)(3p + 2)(5p + 3){7p + 5) 
p(p(p(p(8662p + 12621) + 1175) - 11631) - 5337) + 270 



18(2p + l)(3p + l)(3p + 2){5p + 3){7p + 5) 
128p^{p{5p + 19) + 6) 



9(3p + l)(3p + 2)(5p + 3){7p + 5) 
p(p(p(65 - p(970p + 2159)) + 2997) + 1833) + 90 



J-2I02) 

102) (209) 
i-4i-l 102> 



3(2p + l){3p + l)(3p + 2)(5p + 3)(7p + 5) 
16p(p + l){p{p{123p- 203) + 144) + 180) p(269 - p(97p + 29)) - 15 3 

lj_'3lj_2-tv_l 02) + o 02 ) 

9(2p+ l)(3p+ l){3p + 2){5p + 3){7p + 5) 6(3p + l)(5p + 3){7p + 5) 

32p(p + l)''(p(35p + 52) + 60) 
9(2p + l)(3p + l)(3p + 2)(5p + 3){7p + 5) 
(p + l){p{p{518p - 151) - 921) - 90) 



L~.3L_i\<t>2) 
-t-2-£'-l 102) 



i'Lil02) + . 



9(2p+l)(3p+l)(3p + 2){7p + 5) 
_ 2(p + l)^(p(7p(58p + 1) - 723) - 270) 

9(2p + l){3p + l){3p + 2){5p + 3){7p + 5) 
_l_ 4(p + l)3(p(10p- 17) -30) 

9(2p + l)(3p + l)(3p + 2)(5p + 3){7p + 5) 
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